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BACKLASH IN SECON1 ORDER 
FEEDBACK CONTROL SYSTEMS 

1 . Introduction 

Whenever components are mechanically connected by a linkage 
or a gear train, there exists the basic problem of providing smooth 
operation without looseness in the connecting parts of the mechanical 
system. In general this basic problem is not completely solved, and 
the connecting parts are free to move independently over a limited 
range. This phenomenon is commonly called backlash and its effects 
are usually „ but not always, undesirable. In particular, when back- 
lash occurs in a feedback control system it often causes a small 
amplitude "dither 11 or limit cycle, and when the system reaches static 
steady state without a limit cycle the accuracy of the system is 
affected by the backlash. 

In the usual feedback control system, the backlash occurs in 
two places: in the power drive gear train between motor and load, 

and in the gear train between load and output measuring device. It 
is not difficult to establish the pertinent differential equations for 
either of these cases, nor is there any conceptual difficulty in the solu- 
tion of these equations, but considerable labor is required for the 
solution of even the simplest cases. Furthermore, the solution of 
one or two specific cases does not provide perspective for generaliza- 
tions s and the multiplicity of parameters involved makes detailed 
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analyses unattractive. This report deals with some initial results in 
a long range project to study the backlash phenomena quantitatively, 
with the ultimate goal of providing charts which will be useful in the 
engineering analysis and design of feedback control systems involving 
backlash o 

In general, the simplest linear feedback control system has a 
second order differential equation, and is commonly referred to as a 
second order system. Backlash may and does occur in such systems. 
While it is obvious that most feedback control systems are third order 
or higher, the study of backlash in this report is restricted to second 
order systems only, so as to provide a firm foundation for future 
studies while limiting the variable parameters to a variety which makes 
numerical calculations feasible. In like manner the mechanical para- 
meters are restricted to inertia (mass) and viscous friction; couloumb 
friction, though of some practical importance, is not considered quanti- 
tatively and is permitted in a qualitative sense in only one special case. 
It is recognized that the mechanical inertia and viscous friction may 
be subdivided and appear on both sides of the backlash in varying pro- 
portions, such cases are considered in detail. It is also recognized 
that the gear tooth surfaces must hit each other after the meshing gears 
have gone through the dead space (backlash), and this impact may be 
either plastic or elastic; this effect is also studied in detail. The 
mechanical load on the gear train is considered to be only inertia and 
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viscous friction,, This combination contains only one energy storage 
device (inertia) and is therefore a first order load* In some cases 
the mechanical connection, e.g. shaft, may act as a spring, thus pro- 
viding a second order load which is capable of resonance if excited by 
periodic motion of the proper frequency. This case is not considered 
but is reserved for future study. 

The two basic tools for studying the nonlinear effects of backlash 
are the phase plane-phase space technique and the describing function 
method. The describing function method has been applied to a number 
of backlash studies, and it is intended to devote some future efforts to 
an extension of such studies; however, the investigations reported on 
here are based solely on the phase plane-phase space technique. For 
second order systems the phase space is two dimensional, that is, 
it reduces to a phase plane, and nonlinear as well as discontinuous 
systems can be analyzed using the well 1 n graphical methods of 
integration (method of isoclines, Lienard ! s construction, the Delta 
method) which make the phase plane a powerful tool,, This provides 
a convenient starting point for the analysis, though the graphical 
methods do not apply to higher order systems. 

The phase plane technique with graphical integration was used in 
initial studies recorded here but was applied to only the simph 
cases. The results obtained and the procedures developed provided 
sufficient insight to permit formulation of a rather comprehensive 
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digital computer program to study the effects of many parameters on 
the response of a given system. The majority of the data presented 
here was obtained using a digital computer program; however, any of 
the results so obtained could be duplicated by graphical methods. 

The digital computer was programmed to provide a complete solu- 
tion of the differential equations, and sufficient data was computed to 
provide a complete and very accurate phase trajectory. However, for 
most of this initial phase the question of primary interest was the 
development of criteria for the existence of limit cycles, The computer 
program was therefore provided with a subroutine which printed out 
only the maximum position per cycle. This was done to decrease the 
read out and print out times. By a simple auxiliary command, complete 
data for the phase trajectory was made available; however, the com- 
plete phase trajectory was solved for only a restricted number of cases 
because of the time involved. As a result, the limit cycle studies are 
detailed and fairly complete (more data would improve the results but 
should not introduce any startling features); but many questions con- 
cerning response during the transient period have not been completely 
answered. Much of the remaining work (within the restricted scope 
of this initial phase) is just a matter of data reduction. 
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2 . THE NO LOAD CASE 



For some systems one may neglect inertia and viscous friction 
on the output side of the gear train containing the backlash. This is 
arbitrarily called the No Load Case. Two practical variations are 
common. First, the most commonly encountered, is the case of the 
instrument servo which is driving some simple indicating or measur- 
ing device such as the wiper arm on a potentiometer. In this case 
the actual driven load has a negligibly small inertia and viscous fric- 
tion. There is usually some couloumb friction in this load, but it is 
often too small to affect the differential equations though large enough 
to prevent jitter in the low inertia driven load. The second case of 
interest here is that in some larger systems with driven load of appre- 
ciable inertia 5 the power gear train has negligible backlash, but there 
is a measuring device geared to the load, and there is backlash in the 
gear train between load and measuring device. Note that while these 
two cases arise in systems which are physically different, they are 
mathematically equivalent and thus one analysis suffices. 

The block diagram description of the system is given in Fig, 2-1. 
11 it is assumed that a step displacement input is applied to the system 
when the gears are in contact and the system is driving linearly, then 
the basic linear differential equations are. 
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The dynamic response of the system is shown on the phase portrait 

of Fig. 2-2, and equations 2-1 and 2-2 describe the first portion of the 

trajectory from M to N. At point N, E = O and © = 0 ? the motor re- 

verses and runs in the opposite direction, increasing the velocity and 

displacement of the motor shaft. However the output shaft cannot move 

and the error, E., cannot change until the backlash is taken up. The 

operation during this interval is: 

J M 0 M +f M @ M = KE N (2-3) 



from which 




where 



•f= f/ 2 V KJ M andW " 




( 2-4) 

(- l+2fW r t+e' 2 ^ W i' 1r ) 



Also 

’ _ -2fWr»t 

0 M--"2^ (1 - e ) 



(2-5) 



It is readily seen that the horizontal or E-axis is a locus of points 
at which the describing differential equation changes discontinuously , 
and it may then be called a "dividing line." This particular dividing 
line is a "separation line" i.e. as the state point crosses this line the 
physical system goes into the backlash region and the drive member 
is mechanically separated from the load member During the time 
interval when equations 2-4 and 2-5 apply, and the backlash is being 
"taken-up"; the load member does not move, E and E do not change, 
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FIG. 2-2 PHASE TRAJECTORIES OF AN INSTRUMENT SERVO- 
MECHANISM WITH BACKLASH 
Di -DIVIDING LINE FOR RECOMBINATION 
D 2 -DIVIDING LINE FOR SEPARATION 



but © and © do change: On the phase plane of Fig 2-2 the vertical 

mm ° 



line from N to P indicates the first interval of this type. It can be 
argued whether this line is or is not validly a segment of the phase tra- 
jectory, but the physical interpretation is quite clear: E and E are 

constant at point N until the backlash is completely taken up, then as 
the gear teeth hit (with inelastic impact) the load instanteously acquires 
a velocity and the state point jumps discontmuously from N to P Once 
the backlash is taken up the system once more operates linearly and 
the segment of phase trajectory from P to Q is again governed by 
equations 2-1 and 2-2. The point P is one point on another dividing 
line, this is the line at which the backlash is taken up and the system 
physically recombines, i.e. , the mechanical connection is restored. 
Thus the line may be called a backlash n recombination line 0 n This 
line can be constructed accurately only by a point-by-point calculation. 
However, equation 2-5 can be rewritten: 



where t^ is the time interval required to take up the backlash For 



values of the backlash,^ , the recombination line may be represented 
on the phase plane as shown in Fig. 2-3. Only the tangent line is 
accurate (unless detailed computations are used), but knowledge of 
this tangent line permits derivation of a stability criterion. 




N 



( 2 - 6 ) 
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FIG. 2-3 FAMILY OF BACKLASH DIVIDING LINES WITH 
COMMON TANGENT AT THE ORIGIN. 

.4 MAX* TAN' 1 — 1/2 j 
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2 . 1 A Stability Criterion for the No Load Case 



If equations 2-1 and 2-2 are combined and nondimensionalized 
forming a single equation with error, E, as the independent variable, 
this equation can be manipulated and integrated to define the linear 
phase trajectory on the phase plane as* 



a 2 # 2 2 

E +2f w n EE 4* UJ f\ E 



2iT 



-1 



tan 



2 V 1— ? 

= D e 



E+f O n E 
WnE-v/ 1-f * 



Using the linear transformation 



(2-7) 



Y=E+£w n E,X= w n E Vl— $ 



Z Z 2 \J 1-t 

X +Y =D e v 



2JT . -1 7 

tan — 



(2-8) 



(2-9) 



-1 



Letting Y= p sin yp , X= ^ cos yp , \p = tan Y/X 

Tv* 






(2-10) 



p = De 

The backlash dividing lines are located on the X - Y plane by noting 



that 



Y 

X 



O 

_ E+jTa» n E 
U p E^l- f 2 



( 2 - 11 ) 



but the dividing line for separation is the E-axis, for which E=O s and 
for this line 



Y_ 

V 
u. V 



r 



v 



i 



J* 



2 



tan yj, 



( 2 - 12 ) 
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The backlash dividing line for recombination is not readily determined, 



but the tangent to this curve A on Fig. 2-3, is easily transferred: 

max 



Y_ E/ uy 
X E 



(2-13) 



s/i-r 2 V i- $ 



and for E / E = — 1 / 2 

max ° 



so 



Y 2^-1 



X 



2 f V 1— 



= tan 



(2-14) 

max 



which may be expressed as an acute angle 



tan 






max 



1 — 2f 2 

2|\/l_ i- 2 



(2-15) 



In Fig. 2-4 the backlash line for recombination is sketched, not calcu- 
lated* but the basic condition for increasing oscillations is 

Pl" 1 'P 1 = P 2 cos N't , 2 - 16 ) 

Assuming that the angles are measured from the negative real axis in 



Fig 2-4 and that p 


and are known: 




£ $ 1 






D=p i e \/ 1 — i* 2 




(2-17) 




J*(*P2- 1 ) 




„ Vi— J 2 


Vi— ? 2 


(2-18) 



p 2 =De = pi* 



N>2= 7t - N> t 



(2-19) 
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Substituting , 






( 2 - 20 ) 



cosl|^=- cos 



V,.f 2 






V 



/ 1 



( 2 - 21 ) 



The angle may be chosen arbitrarily, and by inspection of Fig. 2-4 
it is apparent that if the required inequality exists it must be most pro- 
nounced when the point (p^ , 1)^) is l° cate d near the origin on the tangent 

line, in which case = U> Substituting this value for UJ * 

1 1 max 1 1 



( 2 - 22 ) 



l£“ e 



f -1 1 - 2 i * 2 -1 £ 

• — tan — - - - 1 — — tan 



Vi-/ 



2 f Vl- t 



\/\ — t 



Equation 2-22 has been plotted in Fig. 2-5, and the stability point is 
marked. This shows that the system with backlash, but no load, will have 
a limit cycle if £ 0.29, and the existence of the limit cycle is inde- 

pendent of the amount of backlash, so long as some backlash is present. 
Obviously the amplitude of the limit cycle will be determined by the amount 
of backlash. 
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3. PHASE PLANE ANALYSIS C F BACKLASH IN A SECOND ORDER 
SYSTEM WITH DIVIDED LOAD AND PLASTIC IMPACT. 

In the more general cases of backlash, the inertia and viscous fric- 
tion on the output side of the gear train are not negligible . For second 
order systems, such cases can be analyzed using phase plane techniques. 
This is true for any subdivision of inertia and friction between drive 
member and driven member, for any gear ratio, and for plastic or 
elastic impact at the gear tooth surface. While actual graphical integra- 
tion to construct the phase trajectories is too laborious to be a tool 
for generalized analysis, understanding the techniques, and actual analysis 
of one or two specific cases helps build appreciation for the nature of the 
problem Graphical solution points out the physical phenomena and 
boundary conditions, and leads naturally into the preparation of digital 
computer programs. 

When the backlash is taken up and the system is operating linearly, 
the equations of motion are exactly the same as for the no load case. 
Considering here that the motor used is specifically a D. C< armature 
driven motor, and using a detailed notation, the equations of motion are: 



( N 2 j + J T ) © +— (N 2 f + f_ ) © = T 

N ' m L c N m L, c 



T= K 1 I 

o 

V— R I + K e 

2 m 

V K 2 

i = (TT — N e c> 



(3-1) 

( 3=2) 
(3-3) 

(3-4) 
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(3-5) 



V = K ( © - © ) 
3 K c 



b- ( N 2 J + J T ) © +4? (N 2 £ + f_ + 

N m L c N m L 



•S v - 2 



N 



R 



) © 



K K © 

+ 1 3 c 

R 

Where 



K 1 K 3 ®R 
R 



( 3 - 6 ) 



K= Motor torque constant 

1= Motor generated voltage constant 
Error measurement constant 
■= K j K 3 /R 
R=armature resistance 
N=gear ratio 

J= inertia of system without load 
J 1 ? load inertia 

£= viscous friction (mechanical) of system without load 
f= viscous friction of load 

F= N f + K_ K N/R = viscous damping of system without 
m m 12/ , , ® 

load but including motor back EMF 

5: £ l /n 

^commanded displacement 
<3= displacement of motor 
^^displacement of load 
©^displacement of combined system. 

To permit convenient numerical manipulation, let N= 1 and 0 = 1, 

R 

then equation 3-6 becomes: 



• * p* p* 

© + ( m + L ) © = ( 

C c 

J + J T 
m L 



K ) (1 - © ) 

c 



(3-7) 



J + J T 
m L 



Equation 3-7 may be manipulated into the isocline equation in the usual 



fashion: 



© = 1 - © 
c c 



N ,+ 



F F t 
m + L 



1 J + J 
m L 



K/ ( J + J T ) 
/ m L 



(3-8) 



where N is defined to be d © /d © and is the slope of the phase 
1 c/ c 

trajectory as it crosses the designated isocline. 
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Equation 3-8 is the normal isocline equation for a linear second order 



feedback control system. It is expressed in terms of motor inertia 
and friction and the load inertia and friction because this subdivision 
becomes important when backlash is considered 

When there is backlash between motor and load, and when this back- 
lash is inside* the closed loop, the system enters the backlash region 
(i.e. , the gears lose contact) when the slope of the load deceleration 
trajectory is exactly equal to that of the linear system. In terms of 
equations, when the load is drifting free- 



J 



L 





(3-9) 



from which 



d© 



L 





(3-10) 



When the system is operating linearly as described by equation 3-8, and 

the state point reaches the isocline which is numerically equal to 

then the gears lose contact, l e. , the backlash separation line is that 

isocline of the linear system for which N = N = — F I J 

1 L JLf / E 

When the backlash separation line is crossed the load drifts freely 
according to the relationship of equation 3-10, which guarantees a phase 
trajectory that is a straight line The motor shaft does not remain 
stationary It moves, but is driven open loop rather than closed loop, 



*It is possible to have backlash in the gear train between motor and 
load but so located that it is between the measurement pick off and the 
load. Then there is no backlash in the closed loop, which operates 
linearly. This case cannot have a limit cycle and is not considered here. 
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with a time varying input. The equations for the motor shaft motion is 
• « « 

J © + F © = K(© o — ©.) (3-11) 

m m m m R L 

The phase plane representation of these equations is shown in Fig, 3-1. 
The separation dividing line is the isocline passing through 21. Equation 
3-11 cannot be represented by an isocline relationship, so first the 
analysis requires a time solution of equations 3-9 and 3-11, The back- 
lash is completely taken up in the reverse direction when the difference 
between the displacement of the load and the displacement of the motor 
is exactlvA, the magnitude of the backlash. In the solution of equations 
3-9 and 3-11 the initial conditions are determined at the separation dividing 
line and are as shown on Fig. 3-2. Transforming equation 3-9 

= O (3-12) 



2 

s © (s) — s © — © 


+ F 


s e (s) — 0 , 


L L 1 L lJ 


L 


_ ]} 





L 

from which 






©. 



h 



®L (S)= Ms + F l /J l ) + . 

Taking the inverse transform 



(3-13) 



F l/ 



e (t) = e + 

L h i f l j l 

Equation 3-11 transforms to 



(3-14) 



J 


© (s) — s © — © 


+ F 


s © ( s) — © 


m 


m m ^ m^ 


m 


m m ^ 







= K 



© (s) 

s L 



(3-15) 
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INERTIA AND FRICTION ON BOTH SIDES OF BACKLASH 



N = - f/J = SEPARATION ISOCLINE; A= TRAJECTORY 
FOR COMBINED SYSTEM; B= TRAJECTORY FOR 
DRIFTING LOAD; C= TRAJECTORY FOR MOTOR 
AFTER SEPARATION. 
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1G. 3-2 RELATIONSHIP BETWEEN © u , 6 MI , © L| and 9 M| 
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since © (t) = u (t) and © (s) = l/s. 

£\ K ^ 



R 



Substituting equations 3-13 for © (s), and noting equivalent initial 

conditons from Fig. 3-2, equation 3-15 may be manipulated to 



m 



© (s) - 

m 



© K{ 1 - © ) 

m j m i 

h s( s +F |J ) + 1 s^(s+ F |J ) 

mm 7 m ml m 



K © 



m . 



J s (s + F t J ) (s + F I J ) 
m L L m 



) m 

Taking the inverse transform of equation 3-1 6 gives: 

J © Ft K (1 - © ) J 

mm — 

© (t) = © + — ( 1 - e J 

m m F 

1 m 



m 



m 



)- 



m ' m 



"F T 
m 



( 3- 1 6) 



( 1 



F t 
m 



m 



F F 
m L 



- e 



F t 
m 



m 



J F - J F 
L m m L 



* ( 



m 



K © 



m 



) - 



F F 
L m 



F t 
m 



m 



r- J J 
m _ L 

t _ — -p— 

m L 



F 1? 



(3-17) 



) 



m L L 

If equations 3-14 and 3-17 are plotted with displacement as ordinate 



and time as abscissa they may be used as a nomograph to evaluate points 
on the backlash dividing line for recontacting of the gear teeth. If both 
are plotted on the same coordinate axes, then at any instant the difference 



between the ordinate values is the relative displacement © - © , 

Lj m 

and whenever © - © = j£\ the backlash is taken up Since this point 

L m ^ 

also defines values for © and for t, differentiation of the proper 

L 

equations permits evaluation of © and © at that instant. 

m 1 j 
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For repeated cycles these time curves need not be replotted, but may 



be used as a nomograph by simply adjusting the ordinate scales. This 
is due to the fact that the only change in the physical situation is a change 
in the initial conditions, and from Fig. 3-2 there is a fixed relationship 
such that i q 1 — © 

tan t|> = = --g (3-18) 



m l 



1 



where tan 4* is known number. Then 1 -© = © tan 4* which (3-19) 

m j m j 

indicates the required adjustment in ordinate scales. A pair of© and 

Li 

© curves is shown on Fig. 3-3. 



m 

At the instant backlash is taken up both motor and load have finite 
velocities . Therefore both the law of conservation of momentum and 
the law of conservation of energy must be considered If the impact is 
defined to be inelastic (plastic) this is equivalent to stating that sufficient 
deformation occurs to dissipate some energy as heat; then the law of 
conservation of energy need not apply to the dynamics of the mechanical 
system, but the law of conservation of momentum still restricts the 
boundary conditions. It follows that after recontact of the gear teeth, 
the requirement that momentum must be conserved is 



o o 

j © T + j © 

L L mm 



( J + 
m 



V 



(3-20) 



where is the instantaneous velocity of the recombined system. Note 
that the points at which the gear teeth recontact define a dividing line on 
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i. 3-31 DISPLACEMENT CHARACTERISITICS OF MOTOR 
D OF LOAD DURING BACKLASH TAKE-UP. 



the phase plane, and the result of the applying the momentum restoration 
defines another dividing line on the phase plane. These lines may be 
established point by point in evaluating the phase trajectory, but no simple 
equation for these lines has been obtained. 

A phase plane plot for a particular case is given on Fig 3-4 This 
plot shows the three dividing lines, and the nomograph of Fig. 3-3 was 
used in computing Fig. 3-4. Note that the system goes into a limit cycle 
along path ABCD, 
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FIGURE 3—4 PHASE TRAJECTORY and LIMIT CYCLE, SERVO WITH BACKLASH 

Combined System 0c + 0.20c + 0c = 0r 

System with load Separated 0.5©m + 0.2©m s 0r-0l 
Load Only 0.5 ©l * 0 

Backlash = 0.3 RAD 



4. COMPUTATION BY DIGITAL COMPUTER 



4. 1 Introduction 

Use of the phase plane techniques of section 3 is instructive but 
prohibitively laborious if a large number of cases is to be investigated. 
Alternatives are analog simulation or digital computation. The indivi- 
dual equations are readily simulated on the analog computer, but instru- 
mentation of the boundary conditions can be difficult. An analysis using 
analog simulation has been carried out elsewhere^ but was restricted 
to the inelastic impact case. It was decided to use a digital computer 
to calculate the desired data, and the computer was to be programmed 
using phase plane concepts so that all of the data required for phase 
plane plots would be available when needed. The first program* was 
prepared for the inelastic impact case only. It operated successfully 
as verified by point checks calculated long hand on the phase plane . 

While the results were accurate, only a limited amount of data was 
computed due to lack of available time . Upon resumption of the pro- 
ject it was decided to reprogram the computer to handle the general 
case of elastic * contact of the gear teeth, including the case of 
plastic impact as a limiting case This program was completed 
and operated satisfactorily. It is given in detail in the appendices. 

The results reported in the remainder of this text is based on com- 
putations made with this program. 

♦Knoll & Narud 
*New 

^Anderson & Luckett 
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TABLE 4=1 SYMBOLS* 



®R 


(ONE)** 


Angular step input (radians). 
Desired load displacement at steady 
state 


© 

c 


(THETAL) 


Output displacement of load (radians). 


0 

© 

c 


( THE TAD L) 


Output velocity of load (radians) 

sec 


E 




Error in radians. E=©_ — © 

R c 


K 


(KMOT CONST) 


Constant of proportionality relating 
position error to torque developed by 
the motor . 


T 

m 




Developed torque of motor. 


t l 




Torque of load. 


T Lm 




Torque of load referred to motor. 


P 


(RHO) 


Number of teeth on gear 1 
Number of teeth on gear 2 


J 

m 


(MOINERT) 


Inertia of motor measured at motor 


] L 


( J LOADVAL) 


Inertia of load measured at load. 


J 




Total inertia measured at motor 
J = Jm + J 

L 


F 

m 


(FMOTOR) 


Friction of motor measured at motor. 


f l 


(FLOAD) 


Friction of load measured at load. 


A 


(DELTA) 


Backlash of gearbox measured in radians 
at the output 


N 

s 




Slope of phase trajectory (combined 
system. 

d© © 

XT A c c 

N S = d©" = S~ 
c c 


*In this 


table some of the symbols are redefined and have different 


meanings than in preceding sections. 


**( ) Te 


rms indicate computer 


mnemonics referred to in Sec. 3, 



Computer Program Development. 
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n l 




Slope of phase trajectory (load floating 
free). 

d© © 

N = c . c 

L d© © 

c c 


.1 ,1 


© , © 




Velocities of motor and load after impact 


me 


e 


(RESTITUT) 


Coefficient of restitution. 


r 


(ZETA) 


System damping coefficient. 




(OMEGAN) 


System natural frequency. 


Wn 


(OMEGANSQ) 


System natural frequency, squared. 




(FLFTPRIN) 




TT 


Load riction ratio with respect 


* T 




to friction of the system. 


J 

m 


(INERTRAT) 


Inertia ratio. 


J L 


© 

m 


(THETAM) 


Displacement of motor (radians) . 


© 

m 


(THETADM) 


Velocity of motor {radians] 

V sec / 
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4 c 2 Development of Equations 



The equations developed were those of a phase plane analysis. 

Prior to defining the net system equations , the assumptions under which 
the analysis was made will be stated. When necessary these assump- 
tions will be amplified and referred to later in this work 
It was assumed that: 

1 . The gear teeth were initially in contact and the initial 

conditions of the system were all equal to zero. This was 
later proved to be an unnecessary limitation for the study 
of steady state response. 

2. Plastic deformation of the gear teeth during steady contact 
and impact and any torsional deformations of driving shafts 
are negligible . 

3. The inertias of the gears and drive shafts are considered 
as part of the load or motor inertia depending on their 
attachment in the system. 

4. The law of conservation of energy was completely satisfied 
in only the perfect elastic case by maintaining the total 
mechanical rotational energy of the system constant at the 
instant prior to and after impact. The law of conservation 

of momentum is satisfied in plastic, elastic, and intermediate 
cases. When the law of conservation of momentum is the only 
equation required to be satisfied, the energy lost from the 
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system is dissipated in the heat of infinitesimal deformations 



of the gear teeth 

The gear teeth are in contact only instantaneously during 
impact for the elastic and intermediate cases (excluding 
plastic) and that the impulse torques of drive, friction and 
bearing supports, etc. , are zero during impact, e.g. 

5 



T dt = O and © = 0 
m 



dt o 



6. The coefficient of restitution of the two opposing gear 

teeth is the same or is described by an equivalent coefficient 
if the two gear teeth are of unequal coefficients. 

7. Backlash is assumed to be equal at all points on the gear 
circumference. Backlash is measured at the output shaft 

A block diagram of the system considered is presented in Fig. 4-1, 
and symbols have been defined in Table 4-1. 

The equations for the motor and load in contact are* 



<D 

1 

1 

® 

II 

H 




(4-1) 


T = 


R c 

KE = K (e D — 


© ) 


(4-2) 


m 


R 


c 


T = 


O G O 

J T © + f T © 




(4-3) 


L 


L c L c 




T Lm 


-? T l=P < J L ©c tf L®c> 


(4-4) 


© = 


0 © forward 


© = o © + A backward 


(4-5) 


c 


r m 

drive 


c t m drive 
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Fig. 4-1 Block diagram of second order servo- 
mechanism with backlash in the gears between 
motor and load. 
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(4-6) 



0 = p 0- forward and backward drive 

c r m 

• « o 

hence T = p (J © + f_ p © ) 

Lm r Lp m L, “ m 

=p 2 (J. © + f © ) 
r L m L/ m 



T =J ©+£ © +T 

m m m m m Lm 



J © f © p (J © f T © ) 

m c_ + rn c_ + 1 L c_ + L c_ 



P 



(J 



m 



pJ T )e (f pf ) © 

' L r m * 1 j c 



L c + _m + 

P 



T = K (© — © ) 

m R c 



K© =( — +pJ T ) © +(— - + pf T ) © + K© 

R'p r L c p v L c c 

f 



© + 






C (J m P J L) 

+ 

9 



K © 



© + 



c ( !hl + PJ l> 



k© t 



J m P J : 



© + 2 T OJ © + L0 n © =U>n ©T 
c v n c " c i 



where u>n = pK 



J m + P J L 



2 £<*>„ = 



m + 



p 2 l 



J m + P J : 



(4-7) 

(4-8) 

(4-9) 

(4-10 

(4-11) 

(4-12) 

(4-13) 

(4-14) 
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The equation for load alone with gears not in contact is' 



J © + f T © = O 

L c L c 



(4-15) 



or 



© + 



c J . 



© = O 

c 



The equation for the motor alone with the gears not in contact is. 

J © + f © = K ( ©_ — © ) (4-16) 

mm mm R c 

The method of analysis used by New, Ref. e, is used in< order 
to be able to examine the system independently of its Ulp 
Defining first U)^t = t# 



and differentiating Lc/r)dt = dt* 



. , dt* , , , 2 . dt* . 

U) n = ~ — and = ( —JZ ) 



dt 

d© 



© = 



d© dt* 
c 



dt 
d© 



c dt 



dt* dt ’ dt* 



c ** 
= © 



* 



hence © =© UJ 



n 



similarly 



d 2 © 



© = 

c 



dt 



d 2 © 



dt *1 



dt 



2 

d © 



UV 






d 2 © 



(at*) 

finally 



• •?]> «• * 2 

- = © and © = © UJ i 

2 c c c 



'n 



♦ •* dt* 

9 = © U>n = © “77 — 

c c c dt 



Je c d‘=j 



© dt or © = © 

c c c 



(4-17) 

(4-18) 

(4-19) 

(4-20) 

(4-21) 

(4-22) 
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Using the equations of the combined system and the load alone, 

w 0 2 2 

© c + 2 y 0J n © c + ©^ = vjU n ® r combined system (4-23) 

£ l - 

© + — T © = O Load alone (4-24) 

c u c 

L 



and making the indicated substitution to a transformed (*) coordinate 
system for the system equation 

UJn e c +2fu; n © c +u;n© c =Ui n e R (4-25) 

®c + 2 -^ e c + ®c =e R < 4 ' 26 > 

and introducting the slope equations of the phase plane 



N © + 2 {* © + © = © T 

g r ^ r r- T 



c c c R 

If the load equation is first put in the form for the phase plane 



(4-27) 



N © + qL © = O 

L c J c 
J— i 

and then transformed to the (*) coordiante system where 



(4-28) 

0^0 



L 

N © + — — © = O (4-29) 

L c J c 

L 

The same result may be obtained by setting ^n= 1. The results of 
such a transformation require inverse scaling for practical applica- 
tion, examples of which given later. It is noted that the slope of the 
load-free equation is mathematically the same whether in the phase 
plane or in the transformed phase plane. It is pointed out at this 
time that the equations later developed to satisfy the laws of con- 
servation of momentum and energy are independent of the system 
natural frequency. 
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To establi sh equations for the law of the conservation of momentum 

and a relationship satisfying the law of the conservation of energy, 

the schematic of Fig 4-2 is used. It is seen that 

* • 1 

© c = — p © j > ©^- — p ©i and since p = 1 an equivalent 

schematic is shown in Fig. 4-2b. By representing the inertia of the 
load and motor as inertia of the gears figure 4-2c is obtained at the 
instant of impact. The torque equation may be expressed 



T= J © 



= J 



d© 

dt 



Impulse is equal to the time rate of change of momentum. The ex- 
pressions for the rates of change of momentum of the gears treated 



separately may be written: 



-i'll, 'A 

- <s 



r j dt = 



© 

r m j 

)e m 



d © 

IP dt 



dt 



1 J 



dt = J (© _ © ) 

mm m 



dt-^o 



(4-30) 



(4-31) 



and 






dt -40 



dt = J 



m 



, * I 

( © — © ) 
m m 



5 



2 J 



F 

dt -4 o 



r 2 dt = 



for the motor 
d© 

2 1 

p J , dt 



L dt 



(4-32) 



(4-33) 



s 



F_. 

2 J 

dt -40 



dt = 



J 

2 L 



I 






(4-34) 
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Fig. 4- 2a Schematic diagram of motor and load 
inertias with interconnecting gears. 




Fig. 4- 2b Schematic diagram of mechanical system 
equivalent to Fig, 4-2a. 
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Fig. 4- 2c Gear teeth and loads involved 
at instant of impact. 
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where the primed velocities are those following impact. 

If the momentum is to be conserved in the system the impulse 
functions are equal and opposite 



m 



r( © - 



m 



© ) 

m 





and 



(4-35) 



m 



1 • ' 
— (0 - 
p m 



©)=(©- 
m c 



• i 

0 ) 
c 



from substitution of the equations 



(4-36) 



r l = P 2 , V ' P ®1 

The assumption that all other impulse functicn s are zero at the instant 
of impact is restated at this point. 

A definition of the coefficient of restitution will be made for pur- 
poses of this work: 



A 
e - 



• i i 

© © , 
m + 1 



( 



* \ ‘i 

p © - © V 

m c / 



© + ©, 
m 1 



p © - © 

m c 



(4-38) 

It will be shown that this definition of e will satisfy the law of 

conservation of energy. Consider the case of e=0, 

O = — ( p © — © ! ) , p© =© (4-39) 

r m c r c 

in which case for plastic impact (e=0) the gears are moving at the 
same velocity following impact. 

The expression for the coefficient of restitution is examined for 
the case of perfect elastic contact, e = l, for which the law of conser- 
vation of energy is satisfied. 
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(4-40) 



Solution of the restitution equation for e=l yields, 
p ©_ — e r = — P 

* JY) r m r> 



(0 + © ) = © + © 

r m m c c 

For conservation of energy 
i - 2 - 2 i -' 2 - ,2 

— ( j © + j e ) = — ( j © + j © > 

c mm_Lc l mm L, c 



or 



(4-41) 



(4-42) 



m 



( e - © ) = (©-©) 

mm c c 



and factoring 
J 



m • • * • • * » i *» • i • 

— (© + © )(© — © ) = (© + ©)(© - ©) 
. mm m m cc cc 



(4-43) 



(4-44) 



The momentum equation which is independent of e is : 



m 



_ (© -©)=(©-©) 
p J mm c c 

-Lj 

or rearranging J 



( © -©) = (©-©) 
p J m m c c 



(4-45) 

(4-46) 



The factored energy equation is now divided by the momentum 
equation to yield 



P (© +©) = (©+©) 

'mm c r 



I 

) 

m m e c 

which checks with the definition of the coefficient of restitution 



(4-47) 



for e=l . 

The momentum equation implied from the impulse approach that 

i i 

dt o hence © = © and © = © • The definition of the coefficient 

c c mm 

of restitution is independent of position. It may be noted that the 
equations for energy and momentum depend only on the inertia dis- 
tribution, gear ratio and angular velocities. 
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In manipulation of these equations it is seen that load inertia 

2 

may be transferred to the , shaft by multiplication p , the same 

* 

transfer may be accomplished with load friction. Since ©^ , and ©^ , 

are related to © and © it may be reasoned that the results for 
J-j J i 

p= 1 may be extrapolated to physical systems where p / 1 

Since a total energy accounting may be made for the system with 
the equations used, the rotational energy lost is attributed to the heat 
of deformation of the gear teeth. This transfer of energy could m 
fact be determined from the equations used. Thus the conservation 
of rotational energy e= 1 is a special case of a broad interpretation 
of the law of conservation of energy. 
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4 . 3 C omputer Program Development 

The physical equations of the net systems were programmed for 
solution using the Control Data Corporation 1604 high speed digital 
computer utilizing paper tape program input and magnetic tape output. 
An IBM 717 line printer was used to extract data from the magnetic 
tape The Control Data Corporation machine library and the U. S. 
Naval Postgraduate School computer subroutine library were used for 
assembly, Runge-Kutta- Gill numerical integration, and decimal 
output 

Several changes were made in the forms of the physical equations 
of the net system in order to eliminate duplicate computing operations 
and to fit the equations to a form suited to the variable parameters. 

In the table of symbols, computer mnemonic (m) terms which are used 
in the assembly subroutine have been indicated by parenthises and will 
be defined when encountered. 

The equation of the system with motor and load combined 




(4-48) 



was put in the form 




where 



(4-49) 



THETAL = 0 



QNE = © 



R 



c 



THETALD = U = © 



c 




UDOT = © 



c 



OMEGANSO =U) r 2 



41 



(RUNGE) , numerical integration method The RUNGE sub- 
routine required the definition of synonymous terms for the four 
iterative cycles used to produce one extrapolated set of variables. 
The increment of the independent time variable chosen was 0.01 sec. 
for all^“^0 . 1 and 0 . 004 sec for all^* ^ 0.1. 



(4-50) 



(4-51) 



The load free equation 

f L • 

was put in he form @ + - — © = O 

c J T c 
-L 

for the computer 0 + B @ _ 0 

c — c 

A separate set of R UNGE synonyms was used for the load free equation 
in order not to destroy the previous computations before they were 
determined of no further value in computing and also in the making of 

program decisions 

* •* 

© was designated V and ©^ was designated VDOT for the computer 
program . 

The equation for the motor alone 
f 



© + — — © + — _ © = JS_ © 

m J m J c J R 

m mm 

was put in the form 
f 

* t 

© = _ 



(4-52) 



m 



m 



e + f < e _ e ) 

J m J R c 

m m 



D © + C ( ©„ 

— m — R 



© ) 

c 



(4-53) 

(4-54) 



and for RUNGE, the terms 



© = THETAM 

m 

and © = WDOT were used, 

m 



© = THE T ADM = W 

m — 
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The equation representing the law of conservation of momentum 



and the definition of coefficient of restitution were combined to the 

2 

T // 

forms 



© = 

c 



m 



V P j l l 



p © ( 1 + e) + © 

m 



‘.CM 

' m ' 



J 



(4-55) 



• I 
0 

m c 



© =© . — ep© + e © 

r m c 



(4-56) 



for the computer program No provision for additional m terms was 

made to denote the primed values. The additional terms: 

2 

e = RESTITUT p = RHO and p ~ RHOSQ were used , 

Two equations were used to define the boundaries of operation of 

the load and motor when they were acting separately, 

© = © and © = p © + A (4-57) 

c » m c I m 

Two major decisions of the computing cycle were the determination 
of the point where the motor and load would float free, F LOATEST , 
and the response of the system to impact of the load and motor gears 
when the boundary conditions on position were met, COM3 TEST 0 

The first decision, the point of float free, is made on the equality 
or inequality of slopes in the phase plane, i.e* , at separation 



N = N Equating the slopes: 

L s 2 

u n 

__(O R ~0 c ) + B_ = O 



(4-58) 



The above equation is solved after computing each point in the 
combined phase trajectory. If N^* O the system remains com- 

bined. If N + N < O the sytem separates, having the response 
S L " Sh 

motor alone and load alone immediately thereafter. 
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The second major decision, that of whether after impact the 
response should be that of the combined system or the response of 
the load and motor acting separately, is made on the basis of the 
resultant velocities after impact and upon which side of the backlash 
the motor and load positions are found. 

The computer program, its flow diagram, a description of the pro- 
gram operation, and an explanation of the readout of information with 
examples are presented in the Appendix, 
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4 . 4 Methods and Scope of Examination 



Two modes of operation of the computer program were used 
The first mode considered was the phase trajectory. This mode was 
used to check the operation of the computing cycle. Time, velocity 
and position of the load were printed when the system was combined 
When the system was separated, the values of motor velocity and motor 
position were also printed. Although the time increment used for com- 
putations was either 0,01 sec. or 0.004 sec. , printed outputs for the 
first mode were taken only every 0. 1 sec. Additional printed outputs 
were taken at the time of contact of the gears just prior to and 
immediately following the solution of the momentum and restititiL;.'* 
equations. A typical printout of the phase trajectory mode of operation 
is shown in the Appendix . 

The second mode of operation was used for priv Ang only the maxi- 
mum computed positive load overshoot position for each cycle, the 
associated load velocity and the exact problem time of the computation. 
A sufficient number of print- outs for each problem was obtained to 
prove limit cycle existence and average size Since this investigation 
was concerned primarily with the steady state response the second 
mode of operation was the one utilized to obtain the majority of the data. 
This read-out method markedly decreased the data reduction time for 
the problem analysis, as opposed to the analysis of steady conditions 
provided by a full phase trajectory print-out Since the time required 
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for computer read-out was several magnitudes greater than the 
computing time, valuable computer time was saved by this mode of 
operation. An example of the printed output obtained when only maxi- 
mum load position overshoot was of interest is given in the Ajoendix. 
Approximately 1400 phase trajectories were solved using this mode of 
operation. Each solution required an average of three minutes of com- 
puter time. 

For the purpose of general examination of the problem it was 

2 



assumed that p = 1, K= l,U) n = 1, p J^= l r Specific solutions 

were made with p ^ 1, 1 to determine scaling effects. It was 

determined that velocities were scaled in direct proportion with uJ^ 

and time is scaled in an inverse ratio with u) It was also determined 
2 f L 2 J L 

that p ^ — and p — — were valid parameters and in agreement with 
* T m 

the non dimensional parameters used for the investigation. 

The program can also be made to solve linear systems by setting 



= O. 



The major variable parameters used in this investigation were: 



0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0 8, 1.0 

2 f 

P T 

■==■= 0 , 0.2, 0.4. 0.6, 0.8, 1.0 

f t 



2 

P 



m 



= 0 . 1 . 

0.9, 



0.2, 0,5, 0,8, 0.9 
0.8, 0.5, 0.2, 0.1 



= 0 . 3 radians 



e =0,0.6,08,10 
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Additional parameters of 



p 2 < 



= 0 . 04, 0,1, 0 . 9; and ^= 0 01,0.03, 0,1,0.15 radians were 



used in certain instances to examine particular characteristics of the 
response . 

The program was limited somewhat in that the value of 




y— =Oo and zero were excluded due to generation of undefined 
m 

mathematical quantities. 

The value of backlash was made abnormally large, ^ =0.3 radians, 
to allow an easier interpretation of the non-linear response. S-ince 
the influence of^was linear, this caused no inaccuracies. All 
graphs for limit cycle size were plotted with ^ = 0,3 radians. The 
effect of values for other than 0, 3 is also shown in graphical form. 
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5. ANALYSIS CF AFSULTS 

The primary purpose of this study was the determination of the 
conditions for existence of a limit cycle. Due to the nature of the 
computations an analysis of the results of each numeric case pro- 
vided an individual n yes M or "no" answer, i e , a limit cycle existed 
or did not exist. This information plots conveniently as individual 
points, and has been recorded on Figs. 5-1 and 5-2 in condensed, 
semi- quantitative form to give a broad view of limit cycle conditions. 
On Figs 5-1 and 5-2 the small circles merely indicate points at 
which a computation was made, the curved lines were averaged in 
as accurately as possible with the available data. In general these 
figures are sufficiently accurate for an engineering check, since 
parameter values which are close to a boundary obviously lead to 
a long- duration oscillation even if there is no limit cycle. 

The data in Figs. 5-1 and 5-2 can be presented in the form of a 
three dimensional picture which provides considerable perspective. 
Such a presentation is given on Fig. 5-3, which is a three dimensional 
sketch of the limit cycle boundary for the particular case of e= 1 . 0 
The surface of this contour is the limit cycle boundary. For com- 
binations of parameters specifying a point on or below the surface 
there is no limit cycle. For combinations of parameters specifying 
a point above the surface a limit cycle is guaranteed, and the 



48 




Ft 4 ** 4-1 
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Limit Cycle. Contour Chdrll 

(not to SCRt-E. SEE PtCr. 5-&) 
CijclC — O radians 



amplitude of the limit cycle is related to the distance from the surface 
to the pertinent point. This is further elaborated on by the sketches 
of Fig. 5-4. 

For calculation purposes detailed and quantitative curves are 
desired, and the three dimensional figure is not a convenient tool. 
Therefore the results of the digital computer study are summarized 
using a series of families of curves. These are given in Figs. 5-5 
through 5-28 to show the effects of parameter variation on both the 
existence and amplitude of the limit cycle. This is essentially the 
data which was summarized in Figs. 5-1 and 5-2 (limit cycle exist- 
ence), but is more complete and more accurate. 

In all of the computations reported to this point the backlash was 
set a^/\ =0. 3 radians. It was tacitly assumed that all amplitude pheno- 
mena would be directly proportional to the amount of backlash. Since 
this condition was easily checked by additional computations it was 

decided to verify the assumption. The results are shown on Fig. 5-29 

# 

and 5 30, from which it is evident that the amplitude of the limit cycle 
is indeed directly proportional to the amount of backlash. 

Included for purposes of general information are two phase tra- 
jectory plots shown on Fig. 5-31 and 5-32. These show the type of 
information available from the digital computer data. 
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Limit Cycle, ( rad.) 





Conata-nt Profile 

'tS - 0-1' 



F ‘ 3 * 5 “4 

Limit Cycle Contour Chart 
Profile Views 
e = to 

a = 0.3 
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Limit Cycle Error (Radians) 




Fig. 5-5 

Maximum Error of Limit Cycle from Unit Step Input 
Limit Cycle Error (Radians) ( \ for /— — - Variable 

\ A ; Jn i 
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Ft 



4 ??- + 1 

/**fu 



= O 



= O 
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Limit Cycle Error (Radians) 




Fig. 5-6 

Maximum Error of Limit Cycle from Unit Step Input 

a I 

Variable 

Jm 



Limit Cycle Error (Radians) 



0-3 

A 



for 



/2_£l „ 

Ft 



£m_ + i 

/=*fu 1 



e = O 



No limit cycle exists for 



= 0.2 



^ 0 ,£> 



Limit Cycle Error (Radians) 




Fig. 5-7 

Maximum Error of Limit Cycle from Unit Step Input 
Limit Cycle Error (Radians ) ( Q - : ^ — \ for ^ — Variable 

\ A ; Jm 



rr 



fm 
/•*fu 

e = O 



+ 1 



« 0.4- 



No limit cycle exists for Jf 3 ^ 0-5 
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Limit Cycle Error (Radians) 




Fig. 5-3 



Maximum Error of Limit Cycle from Unit Step Input 
Limit Cycle Error (Radians) ^ 0^ — 'j 

/*h. - t±- 



Jm 



Variable 



rr 



/»*f L 
e = O 



+ 1 



= 0.6 



No limit cycle exists for 
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0.5 



Limit Cycle Error (Radians) 




Fig. 5’ 9 

Maximum Error of Limit Cycle from Unit Step Input 



Limit Cycle Error (Radians ) f Or— \ for Z '— , Variable 
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A . 



rr 
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= 0.6 



= o 



No limit cycle exists for ^ 0-6 
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Maximum Error of Limit Cycle from Unit Step Input 
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l A / J-tn 
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/Lh. = 

Ft 
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No limit cycle exists for Jf 5 
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= 1.0 
^ 1.0 



Limit Cycle Error (Radians) 










/ziii 

J ryt 



Fig. 6"il 

Maximum Error of Limit Cycle from Unit Step Input 



Limit Cycle Error (Radians) ( \ for Variable 

\ A / Jm 



/£l± = 

Ft 



1 



£ 

- 3 * 

= 0.6 



= O 
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Limit Cycle Error (Radians) 




Fig. 5-lZ 

Maximum Error of Limit Cycle from Unit Step Input 



Limit Cycle Error (Radians) ( for /— Variable 

\ A ) Jtn 
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Ft 



/Lh. , 



1 
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e =0.6) 
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Maximum Error of Limit Cycle from Unit Step Input 
Limit Cycle Error (Radians) f ^^-Wor — - Variable 

\ A J Jrr. 



rr 



f"i- + 1 
/ 3 *‘fL 1 



= 0.4- 



e = 0.6 

No limit cycle exists for 0.5 
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Limit Cycle Error (Radians) 
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Limit Cycle Error (Badlans) 



for 



Maximum Error of Limit Cycle from Unit Step Input 

Jm 
Ft 






Variable 



±"*_ jg. i 

l 1 

= o.a 



= 0.6 



No limit cycle exists for Jf 5 ^ 0. 5 
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Limit Cycle Error (Radians) 




Fig. 5-15 

Maximum Error of Limit Cycle from Unit Step Input 
Limit Cycle Error (Radians) for Variable 

. , 

Ft 



e = 0.6 



+ 1 



= 0.3 



No limit cycle exists for y* ^ O-Q 









Limit Cycle Error (Radians) 
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Maximum Error of Limit Cycle from Unit Step Input 

Limit Cycle Error (Radians) ( ^ for S'— Variable 

\ A ; jtn. 



/!!£l , 

Ft 



+ 1 

/’'fu 
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e = O. 6 
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Limit Cycle Error (Radians) 




Fig. 5-17 

Maximum Error of Limit Cycle from Unit Step Input 
Limit Cycle Error (Radians) ( for Variable 

\ A J Jw 




o 






0.6 



Limit Cycle Error (Radians) 




Fig. 5-18 

Maximum Error of Limit Cycle from Unit Step Input 
Limit Cycle Error (Radians ) ( for / ' — Variable 

\ A | Jm 

/tK . ■ 



rr 



+ 1 



= 0.8 



/»*f L 
e =0.8 

No l imi t cycle exists for y* ^ O.Q 







Limit Cycle Error (Radians) 
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Maximum Error of Limit Cycle from Unit Step Input 
Limit Cycle Error (Badians) ( — for ^ - L - Variable 

1 A j Jm 

- 



Ft 



/ 5 *fL 

e = O.Q 



+ 1 



= 0.4- 



No limit cycle exists for Jf* ^ 0.5 









Limit Cycle Error (Radians) 




Fig. 5-£0 

Maximum Error of Limit Cycle from Unit Step Input 

/£A 



Limit Cycle Error (Radians) ( ■ ■ j for — - L 

Jtn 



/Z^L - 

' Ft 



A 

i 



Variable 



la. + 1 
/••ft 1 

= 0.3 



= 0.6 



No limit cycle exists for Jf 3 ^ 0.5 
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Limit Cycle Error (Radians) 




Fig. 5-ZI 



Maximum Error of Limit Cycle from Unit Step Input 

0 



Limit Cycle Error (Badians) ( \ for — Variable 

A j Jm 



/££l 

Ft 



fnn. 

/»*fu 



+ 1 



= O.Q 



e = 0-6 

No limit cycle exists for ^ 0.6 
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Limit Cycle Error (Radians) 




Fig. 

Maximum Error of Limit Cycle from Unit Step Input 
Limit Cycle Error (Radians) ( for / — — Variable 

\ A / Jm 



rr 



fra. +. t 



e - 0.8 

No limit cycle exists for 1-0 



$ 



i 



Limit Cycle Error (Radians) 




Fig. 5"&3 



Maximum Error of Limit Cycle from Unit Step Input 

| 

for /— — — Variable 



Limit Cycle Error (Radians) 

A = 

Ft 



/’•fu 1 



Jtn 

= O 



= 10 
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Limit Cycle Error (Radians) 




Fig. 5-£4 

Maximum Error of Limit Cycle from Unit Step Input 
Limit Cycle Error (Radians) ( for Variable 

\ A J Jtn 

/fih. , - 



rr 



+ 1 



= 0 . 5 . 



e = t.O 

Wo limit cycle exists for J ^ 0.& 



Limit Cycle Error (Radians) 




Fig. 5~^5 

Maximum Error of Limit Cycle from Unit Step Input 
Limit Cycle Error (Radians) ( for / J ' Variable 

\ A ) Jfn 

/£l± . -r±- 

v nrt 



rr 



/°*f I. 

e = 1.0 



= 0.4 



No limit cycle exists for f 2:0.6 
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Limit Cycle Error (Radians;) 




Fig. 5-^6 

V 

Maximum Error of Limit Cycle from Unit Step Input 



Limit Cycle Error (Radians) 
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/L A = 

Ft 






Jtn 



Variable 



£irt 



+ 1 



= 0.6 



e = 1.0 

No limit cycle exists for ^ O- 6 
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Limit Cycle Error (Radians) 
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/ 

Maximum Error of Limit Cycle from Unit Step Input 
Limit Cycle Error (Radians) ( for ■ ' Variable 

\ A ) 



Ft 



A J 

^ 

-r + i 



£m_ 

/’'fL 



= 0.6 



e = 1.0 

No l imi t cycle exists for Jf* O.S 

^6 






/ 



Limit Cycle Error (Radians) 




Fig. 5 -eL& 

/ 

Maximum Error of Limit Cycle from Unit Step Input 

Variable 

'1 A | J*n 

2 



Limit Cycle Error (Radians ) ( ^77-^ for 

\ A ) Jmi 

= f.O 



/LK = 

Ft 



fm 

/»*ft 



+ 1 

e = 1,0 



No limit cycle exists for J> 1.0 
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f 






6. APPLICATION OF RESULTS 

In this section, the results are applied, to hypothetical physical 
system s . Methods are indicated by which physical systems may be 
analyzed or designed. 

It is desirable to find an equation relating certain system variables, 
therefore the equations’ 



UJ n - 



and 



p K p K 

2 

J + p J J 
m r L m 



2 

1 + p J. 






£ + of 

m r 



? £ 



m 



m 

i 






+ 1 



J m + P J L 



^ + P 2 J. 

m L, 



will be mampalated to the form 

f 



m 






2 

P f 



m 



P 2 f 



+ 1 



' /J m s/‘ 



+ P J 



r p 



m 



2 _ - 



1 + p J 



or 



\/p v K 



m 



p 2f i 



m 



+ 1 



J K (l + p 2 J \ 
m ^ r L / 



m 



m 



1/2 



( 6 - 1 ) 



(6-2) 



(6-3) 



(6-4) 



In view of the form of Equation 6-4 for and recalling that the 

r 2 2 

variable parameters used on the figures were p J and p f = 

L L 



m 
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. 






J 



1 



it is seen that the expression for 



2 

P 



£ 



L 



f 

m 



+ P 




f 

m 




+ 1 



in fact includes the other two variable parameters. 

The value of overshoot and hence limit cycle are not affected 
by Ul 0 ; however, the time and velocity are scaled for UJ n / 1.0. Thus 
the transient response characteristics may be related to the same 
parameters used in these applications. 

If a motor and load are selected for a system and the limitation 
is imposed that system operation shall not result in a limit cycle. 
Equation 6-4 and the limit cycle existence charts may be used to 
determine the unspecified system parameters , K and p. The 
coefficient of restitution for the proposed gear train material must 
also be determined. 

To proceed, select a value of K and p, then determine the value 

of £ from Equation 6-4. Enter the limit cycle existence chart for 

. 2 

the value of e, system and the parameters p and 

T~ 



f 

m 

Tt, 



+ 1 



Determine if the selected parameters describe a point of no limit 
cycle. Since limit cycle existence is independent of backlash O 

excluded), backlash is not a variable parameter. Successive trial 
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.. 



values of the unspecified parameters may be required to meet the 
limitation of no limit cycle. If it is also required that certain transient 
response characteristics be met, these conditions cc ild be examined 



by a similar method. 

In the event that no acceptable parameters are found to satisfy 
the system requirements for no limit cycle, the figures for limit cycle 
size may be examined to obtain a minimum size limit cycle. The 
method of solution is similar to that previously described. 

Example: The material intended for use in the gears has an 

e= 0.6. The motor and load have the parameters f = 0 . 64, J =0.25, 



m 



m 



f = 0.64, J =1.0 and it desired that the system operate with no 
L L 

limit cycle. If the values are selected K-4.0, p= 0.5, the values 



p J T =1.0 and f =4.0, 
1 > m 



= 0.2 result. 



m 



m 



p f i 



+ 1 



p f i 



3/2 



Solve Equation 6-4 for = . 5 ( . 64) ( 5+ 1) 

[25 (4) (l+lj 1/2 

= 0.4. Enter the chart for limit cycle existence, Fig. 5-1 
e=0. 6,^* = 0.4, 



, thus 



f 



=0.2, p J = 1.0 



m 



P 2 *, 



+ 1 



m 



It can be seen that limit cycle will exist for these parameters. 
Since K is most easily varied in equation 6-4, a new K is selected 
at 1 .0 in an effort to find parameters for no limit cycle. It is seen 
that 



C = C A = 0.4 (2) 

3 K = 1 .0 3 K=4.0 N 
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with all other parameters constant, thus£* =0.8 Enter Fig. 5-1 
with e= 0 . 6 



£ = 0.8 ^ = 1.0 1 

^ J 



f 



= 0.2 



m 



m 



p 2f i 



+ 1 



These parameters describe system operation with no limit cycle and 
are well within the area limits. The transient response of this system 
may be undesirable and various combinations of p and K would have 
to be tested to satisfy the additional requirements of the problem. It 
will be restated that the existence areas of Figs. 5-1 and 5-2 are not 
exact. 

The figures of limit cycle size may be used in the same manner as 
the existence charts. The same parameters of the first case are 



assumed : k=4.0, p=0.5, pj = 1.0, 

L 



1 



= 0.2 



m 



m 

p 2f 



+ 1 



f=0.4 - -L 

e=0.6,A = 0.1 radians. Enter Fig. 5- 12 with the parameters given. 
The ordinate, Limit Cycle Error ( 0.3 ) =0.13 radians, is obtained. 

A 

The maximum positive error in the limit cycle = 0.13 (0.1 )= 0.0434 

0.3 

radians. Enter the same figure for e=0.6,^*=0.8 



?\ - >-° 



f 



= 0.2, A =0.5 radians; this is the second 



m 



m 



P 2f i 



+ 1 



case 



of examinati on f or existence. No ^ = 0.8 curve exists; hence no 
limit cycle exists for this choice of parameters for any nominal value 
of backlash. 



85 



7 



• > A C; TIC IMPACT 



T r> a -ct^tt nr ^PONsr- CHArx.; - cterictics, ^ 

The effects of backlash on the transient response to a step displace- 
ment input is also of interest, so a relatively exhaustive study was made 
of the plastic impact case. The char cteristics of interest are the peak 
overshoot Mpt, the settling time and the peak time. Of these, the peak 
overshoot, Mpt, produced the most satisfactory pattern with the result 
that d< ta charts and equations have been prepared to assist in analysis 
and design Settling time and peak time were much less satisfacv . y to 
work with The results obtained to date are presented here, but more 
work remains to be done . 

The bulk of the computed data assumes that the backlash is 0. 3 
radians To use this data with other values of backlash it is first necessary 
to eetablish the effect of varying the amount of backlash. This is done on 
Figs. 7-1 *2,3. These figures merely verify the intuitively obvious 
fact that the amount of overshoot must vary linearly with the backlash, 
but that the proportionality constant depends on the various parameter 
values 

7.1 Maximum Overshoot; Plastic Impact Case 

Figs. 7-4, 5, 6,7, 8, and 9 present the results of computer calcu- 
lations in the form of families of curves. Each family is a plot of 
2 2 f * 

Mpt vs p J / J for a specific p /FT ratio with £ as the family 

LI m 



K & A Figs. 5,6,7 

Figs. 10, 11 , 12, 13, 14, 15 in K & A 
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parameter. All are for the case of plastic impact These curves can 
be used to predict the peak overshoot of any arbitary second order system 
with backlash and plastic impact by reading off the pertinent values and 
inserting in the equation 



Mpt 






A = A i Mpt 


A=0 + 0.3 


Mpt 


A= .3 ~ M Pt 




sr s si 




£ ~ S’ \ 



A=D 

S=Si 



(7-1) 



As was pointed out earlier, the solutions were made and the results 

were plotted using a backlash angle of 0 3 radian. Thus the values of 

amplitude of maximum overshoot on the curves are significant only in that 

they represent the slope of the overshoot versus backlash straight lines. 

This slope must be calculated in order to scale the value of overshoot read 

from the curve for the magnitude of backlash actually expected Therefore, 

in order to refine the prediction scheme, other design curves, figures 7-10 

through 7*15, were calculated using the curves of figures 7-4 through 7-9 

dM 

just discussed. Along the curves the slope pt was calculated by, 

d A 



dM 

pt 



dA 




M 

pt 



4=0 






i 



. 3 



(7-2) 



These values of dM were then plotted versus inertia ratios at 

Pt 

constant friction ratios before. These curves give a more straight- 
forward prediction method with the following equation' 



K & A , 16-21 
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(7-3) 



M 



Pt 



= M 



A=A 1 



pt 



A=o 

r=?, 



+ dM / d A| 
pt 



k-c t 



x A 



In order to demonstarte the procedure for using the design charts 
shown in figures 7 4 through 7-15, &n example problem is solved obtaining the 
maximum overshoot, M , for a second order servo with backlash. The 

Pt 

differential equation for the combined motion of an armature controlled 
motor is: 



© + 

c 



f + f 
t 1 



j +p j 

m ' l 



© + 

c 



P K a 



J J T 

m Lj 



p K a 

© = -L-± e (7-4) 

L _ c T . r 



J +P J 
m r L 



Given the following specifications of a D, C „ shunt motor coupled 
through a gear box to a load: 

Horsepower = .008 ? 

Rated RPM = 4000 



R 


282 ohms 


K 


= .233 volts/ radian/ second 


g 


J 


-6 2 
= 8.2 X 10 slug feet 


m 


f 


= 3. 34 X 10‘ 6 lb. -feet/ radian/ 


m 


K 


= 10 volts/ radian 


e 


P 


= i 


A 


= 05 radians 1 


j 


-6 r 2 
= 16.4X10 slug-feet 


L 


f 


= 10.4X10 lb . -feet/ radian/ 


L 


K 


= . 1464 lb. -feet/ampere 
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The required values 



K = 
a 



K K 
t e 

R 



shown in Equation 7-4 are obtained 

■ 1464 x 1 0 _ .0052 LB-FT 

282 



as follows - 




14 5 



RAD 

SEC 



13.74 x 10 



LB-FT 

RAD/SEC 



f t = f T + f i = < 13-74 + 121 ) x 10 



135 x 10" 6 



LB-F T 

RAD/SEC 



K K 

f _ t g _ _ 1464 x .233 _ 121 x 10 

1 R 282 




LB-F T 

RAD/SEC 



2 

P 



J T / J 

1^ m 



(lp 1 6 . 4 x 10 



8.2 x 10 



Tr 



2.0 




10.4 x 10 ~ 6 _ .077 

135 x 10 ' 6 



Two methods are shown for obtaining the maximum overshoot with a 
backlash angle of .05 radians in the gear trains. 

The first method is accomplished using the curves of figures 7-4 
through 7-9 The procedure is simple, but requires more steps than does 
the second. Also one must remember that these curves are drawn for a 
backlash of 0 3 radian and must be scaled for the expected amount of back- 
lash. The second method makes use of the design charts of figures 7-10 
through 7 - 15 . 
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Method 1 



Z Z 

Entering figure 7 4 ( n f / F_ = 0 1) , at o J T /J = 2.0, and a 

1 L / T ( L m 

damping ratio, = 0.19, the amount of overshoot for a backlash angle, = 



0 . 3 rad . is M = 1 . 658 . Then: 

pt 



M 



pt 





= M pt 


+ .05/. 3 


M 

pt 


M 

pt 


<4 = .05 




II 

o 




II 

Od 


f=19 




r= 19 




o 

II 



A=o 

r = ■ 19 



From second order linear characteristics 
M 



hence , 



pt 



pt 



Method 2 



M 



pt 



= 1.545 


=1.545 + 


.05/ . 3 


1.658 - 1 . 545 


= 1 . 545 + 


.0189 




= 1.564 




= M 


+ dM /dA 


^=•05 


pt 


A =0 


r=-i9 




f = .19 



X A 



2 2 

from figure 7 10,(p £ J F T = 0 . 1) , at p J J J m = 2 • 0 - and a dampinj 

ratio J"* =0.19, dM / d/\ = 0.357 

pt 



Also: 



M 



pt 



= 1.545 



A=0 

f =.19 



hence , 



M = 1 . 545 + . 357 X . 05 

pt 

= 1.545+ .018 
= 1 . 563 
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The procedure for using the charts rem? ; ns unchanged for servo 
systems which utilize other types of motors. The only difference which 
arises is in the calculation of the chart parameters from the known 
system constants. For a field control motor if the field time constant 
is negligible, the equation is similar to the armature control motor equation. 



<<» 

© + 

c 



© + 

c 



P K a 



o— © = 

£ c 






© 



J +p J T J +p J T J +p J T 

m » L m l L . m ! L 



R 



( 7 - 5 ) 



In like manner the equation for the system with a two phase motor is 

P 



© + 

c 



F K 
t n 



J + p J 



m 



? J L 



© + 

c 



K pK 

l m 



m 



J m + P J L 



© = 



C J m + P J L 



(7-6) 



where 



K = - 
n 


dT 

dn 


(from 


characteristic 


curves of two phase motor) 


K = ■ 

m 




(from 


characteristic 


curves of two phase motor) 


<3 e 1 



The calculation of the chart parameters is a simple matter once the defin- 



ing equations are known. 
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7 2 Maximum Overshoot, Elas tic Impact Case . 

Investigation of a limited number of cases reveal that the amount of 
overshoot with plastic impact is always greater than that with elastic 
impact. The difference between this maximum overshoot decreases as 



— — increases. However, under no circumstances was the effect 

f t 



onmaximum overshoot greater than .01 radians. The effect on the time 
of maximum overshoot appears negligible. In a few cases the time was 
different for the plastic and the elastic cases but the difference was very 
small and no pattern of difference could be ascertained. 

These results point out the validity of the original assumption made 
for this study. In a practical system there will exist, of course, some 
amount of gear bounce when the gears contact each other. The effect on 
maximum overshoot, however, is so small that the curves obtained m 
this study remain reasonably accurate 
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7.3 Settling Time 



In determining the effect of backlash on the settling time of a second 
order servo system, it was necessary to first find a definition for settling 
time which would be suitable for all of the cases which were investigated. 
Normally, settling time is defined as the time for the output to reach 98% 
of it final value. This definition has real meaning, however, only for 
linear stable systems in which there are no sustained oscillations. For 
systems which do exhibit small sustained oscillations, the settling time is 
sometimes defined as the time required for the amplitude of oscillations to 
decrease to within 2% of the final value. For very small amplitudes of 
oscillations this definition gives an unrealistic value for settling time. 

For the purpose of this study settling time is considered as the time re- 
quired for the system output to settle to 98% of the difference between the 
magnitude of the input and the amplitude of the sustained oscillations. 

For all cases for which the damping coefficient was less than 0.6 
there were no conclusions as to effects of backlash on the system settling 
time since the magnitudes of the sustained oscillations were always quite 
large for the amount of backlash used in the solutions and the settling time 
was difficult to obtain accurately. For those cases for which the damping 
coefficient was greater than 0.6, the sustained oscillations were small 
if they existed and an analysis revealed the following general conclusions: 
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1* As the size of the backlash angle increases, the settling time 
will increase provided the other parameters are kept constant. For the 
cases investigated, this increase in settling time was nearly proportional 
to the backlash angle. 

2 

2 As the ratio n J T /J increases, with other parameters constant, 
* L m r 

the settling time increases. 

2 

3. At the ratio * f / F is increased, with other parameters 
v L t 

constant, the settling time decreases. 

The data included in Tables I - III show these general conclusions 
and also some quantitative effects of backlash on settling time. 
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TABLE I 



SETTLING TIME DATA 

Hr 1 - 0 



f 


“M 


p2f L 

TO ' 


Radians 


Seconds 

*• 


.6 


1/9 


.2 


.3 


28 


.6 


1/9 


.4 


.3 


21 


.6 


1/9 


.6 


.3 


16 


.6 


1/9 


.8 


.3 


13 


.6 


1/9 


1.0 


.3 - 


7 


.6 


1/4 


.2 


.3 


28 


06 


1/4 


.4 


.3 


24 


.6 


1/4 


.6 


.3 


16 


.6 


1/4 


.8 


®3 


14 


.6 


1/4 


1.0 


.3 


9 


.6 


1 


.2 


.3 


29 


.6 


1 


.4 


.3- 


28 


.6 


1 


.6 


.3 


23 


.6 


1 


.8 


.3 


17 


.6 


1 


1.0 


.3 


10 


.6 


4 


.2 


.3 


29 


.6 


4 


.4 


.3 


28 


.6 


4 


.6 


.3 


26 


.6 


4 


.8 


.•3 


23 


.6 


4 


1.0 


.3 


10 


.6 


9 


.2 


.3 


29 


.6 


9 


.4 


.3 


29 


.6 


9 


.6 


.3 


26 


.6 


9 


.8 


.3 


23 


.6 


9 


1.0 


.3 


10 



TABLE II 



SETTLING TIME DATA 



. - 1 .* 0 



f 


/j M 


FT 


Radians 

A 


Seconds 

t 

& 


.6 


1/9 


.8 


.03 


5 


.6 


1/9 


.8 


.10 


7 


.6 


1/9 


.8 


.15 


8 


.6 


1/9 


.8 


.30 


13 


.6 


1 


.8 


.03 


11 


.6 


-i 

JL » 


oQ 


.10 


14 


.6 ' 


n 

± 


.8 


.15 


16 


.6 


1 


.8' 


.30 


17 


.6 


4 


.8 


. .03 


12 


.6 


4 


.8 


.10 


14 


.6 


4 


.8 


.15 


17 


«6 


4 


.8 


.30 


23 


.6 


9 


.8 


..03 


13 


.6 


9 


.8 


.10 


16 


.6 


9 


.8 


.15 


18 


."6 


9 


.8 . 


.30 


23 
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TABLE III 



SETTLING TIME DATA 
LJr» =1*0 



y 


/j K 


FT 


Radians 


Seconds 

t 

B 


.8 


1/9 


.2 


.3 


8 


.8 


1/9 


.4 


.3 


.7 


.8 


1/9 


.6 


.3 


6 


.8 


1/9 


.8 


.3 


5 


.8 


1/9 


1.0 


.3 * 


5 . 


.8 


' 1/4 


.2 


.3 


20 


.8 


1/4 


•4 


.3 


9 


.8 


1/4 


.6 


.3 


7 


.8 


1/4 


.8 


.3 


6 


.8 


1/4 


1.0 


.3 


5 


.8 


1 


.2 


.3 


43 


.8 


1 


.4- 


.3 


31 


.8 


1 


.6 


.3 


10 


.8 


1 


.8 


.3 


8 


.8 


1 


1.0 


.3 


6 


.8 


4 


.2 


.3 


48 


.8 


4 


.4 


.3 


36 


.8 


4 


.6 


.3. 


12 


.8 


4 


.8 


.3 


9 


.8 


4 


1.0 


.3 


6 


.8 


9 


.2 


.3 


32 


.8 


9 


.4 


.3 


43 


.8 . 


9 


.6 


.3 


13 


.8 


9 


.8 


.3 


9 


.8 


9 


1.0 


.3 
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7.4 Time of Maximum Overshoot 



For each of the cases investigated the time of maximum overshoot 
(peak time) was recorded from the printout of the solution. It was 
found that the presence of backlash causes a slight increase in the peak 
time over the linear system for those cases where . 5 , and has 

negligible effect in those where J <^.5. When p^f^ ls increased from 



ze 



ro to 0 5, other parameters constant and^** ^ 5 ^ peak time increases 



slightly and remains essentially constant for p f ^>5. As 



h 

J /J is increased, other parameter s constant and (* . 5 , peak tim< 

T m ^ 



P J L 

is decreased. 
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7.5 Other System Characteristics. 



A careful study of the results as plotted in figures 7-4 through 7-9 
reveals some interesting phenomena. At all friction ratios it is observed 
that as the ratio of load inertia to motor inertia is made smaller the ampli- 
tude of the overshoot is less. There is a point where the overshoot is identi- 
cal to that of a linear second order system, in other words, where backlash 
has no effect on the overshoot Note that this point occurs at different 
inertia ratios for each value of damping ratio. Also this point varies de- 
pending on the value of friction ratio If the inertia ratio is decreased beyond 
this point then the amplitude of overshoot remains constant at the value of 
the linear case Figure 7-16 is a plot of the boundaries of the combinations 
of variables which will cause an effect on the amplitude of overshoot if the 
system has any appreciable amount of backlash The coordinates are fric- 
tion ratio and inertia ratio with the boundaries labeled with values of damping 
ratio. Areas to the right of a boundary are combinations of the variables 
which affect overshoot. These curves were plotted using information deducted 
from figures 7-4 through 7-9. Three points were selected as test points to 
prove the validity of the boundaries as plotted. These points are shown on 
the figures as points n A M , ,! B n , and "C M . The system parameters used at 
the test points were: 

Point A Point B Point C 




F 



T 



= 0.8 






pJjJ =0.25 

p L m 




1 14 



The system was solved on the digital computer to obtain the maximum 
overshoot for these three sets of parameters. According to the curves 
point A falls to the left of the boundary of the r = .8 area. Therefore, 
the overshoe t should be the same as that of the linear case. The com- 
puter solution gave an overshoot of 1.019, which is, in fact, identical 
to that of the linear case. Similarly, point B falls to the right of the 

=0.6 boundary and in the area where the amplitude of overshoot should 
not be the same as for the linear case. The computer solution gave a maxi- 
mum overshoot of 1.110 which is slightly greater than the linear case. 

Point C falls to the right of the j =0.3 boundary and in the area where 
the amplitude of overshoot should be different from that of the linear 
system. The computer solution gave a maximum overshoot of 1.41 while 
that of the linear system is 1.372. From these boundary curves one can- 
not determine the exact effect of backlash on the maximum overshoot of a 
system, but can determine in which systems backlash will cause an appreci- 
able effect. 

Another phenomena which is revealed in the results as plotted in 
figures 7-4 through 7-9 is that for all cases the greatest maximum over- 
shoot occurs in the range 0.1^ p^J ^ 10.0 for the range of 

2 

a J /J tested. However, it appears that since the curves are trending 
f L m 

2 

upwards at the high end of the p J^/J^ axis that at very large inertia 
ratios the overshoot would be greater. It is felt that these large inertia 



ratios are outside the practical range. An attempt was made to explain 



these phenomena but without success. It is often very difficult, if not 
impossible, to find a physical explanation for certain behavior caused by 
non-linearities They can merely be pointed out from analyses of experi- 
mental information and accounted for accordingly 

It is interesting to note that when the ratio of load friction to motor 
friction is below 0.4, as m figures 7-4 and 7-5 and 11, there is a possi- 



dependmg on the inertia ratio. At first glance one would reason that if 
the system is critically damped there should be no overshoot. However, 
it can be seen that if the friction ratio is small and the inertia ratio is 



trajectory, is made small. Therefore, when the load separates from the 
motor, its velocity remains fairly constant at a relatively high value and 
its momentum carries it past the correspondence point. 



bility of an overshoot for the critically 




large then the ratio of F /J , which is the slope of the drifting load 

I-* L 



117 



8. CONCLUSIONS 



Since this report is to be the first of a series, it is apparent that 
the results are incomplete and the conclusions must therefore be brief: 

1 . It has been established by calculation the t a second order feedback 
control system with backlash may or may not exhibit a limit cycle. 

2. The existence of a limit cycle as a function of various parameters 
has been established, and a series of charts suitable for analysis and 
design have been prepared. 

3. The amplitude of the limit cycle depends on the amount of backlash, 
and charts illustrating this dependence have been prepared. 

4. Peak overshoot depends on the system parameters and the amount of 
backlash. Charts and equLtions suitable for analysis and design have 
been prepared 

5. The effect of elastic impact at the gear teeth is significant, but much 
less so than was anticipated. Most of the design charts include this 



effect . 
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STOP 




yes 
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exit' 1 




aTOPAVe 
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APPENDIX A (continued) 



EXPLANATION OF FLOW DIAGRAM BLOCK SYMBOLS 

,1. Clear fault stop. 

2. Print last overshoot stored in LIMITBUF (symbol 888 ). 

Print out position in phase trajectory of .fault stop 
(symbol 666 ). 

5. End of file. 

!<-. Actuate fault stop. 

5 . Load zeta from Index 1 (b^). 

6. Print last overshoot ' stored on LIMITBUF (symbol 888 ) . 

Print position in phase trajectory where it was manually 
stopped (symbol 1961 ) . 

7 . Stop if Zeta = 0. 

8. 2 jr u» . a 

9 . Is zeta = 0.1? 

10. Is zeta less than 0.1? 

11. Set up computation time of 0.01 seconds for problem and 
set up for print at every 0.1 seconds when zeta is greater 
than 0.1. 

12. Set up computation time of 0.004 seconds for problem and set 
up for print at every 0.1 seconds when zeta is equal to or 
less than 0.1. 

13. Set up limit cycle print-out routine for 20 overshoots and 
average the last 8. 

14. Set up number of limit cycle print-outs and averaging routine 
depending upon the size of zeta: 
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f 


Overshoots 


Start 

Averaging 


Divide 
by * 


Comp. 

Time 


^ 0.8 


6 


k 


2.0 


0.01 


0.6 


7 


5. 


2.0 


0.01 


0.5 


8 


5 


5.0 


0.01 


0.4 


10 


7 


3.0 


0.01 


0.3 


14 


10 


4.0 


0.01 


0.2 


20 


12 


8.0 


0.01 


0.1 


17 


12 


5.0 


0.004 



*The last row printed has this number printed on the far left 
if the problem has a normal solution. This number indicates 
how many overshoots were averaged. The average value has 1.0 
subtracted from it and is printed in the far right of the same 
last row. 

15. Page reject. 

16. Clear buffers and set up for new problem. 

17 . Restitution, e, from Index 3,(b,). 

3 

Backlash, ^ , from Index 4 (b^). 

J m from Index 5 (b_ ) . 

5 

Jr * Jm , K Jlir* 

/’* =JLi “ = 0; — 

from Index 6 (bg). 

( f y? T )( ZfU n )(Jr) - fu i = & 

-^7 t -/° 2fL -f" J ^ “ D 

18. . Have all values been cycled through? 

19 . Set STOPNOW after this run to stop absolutely. 



124 



20 



= 0 ? 



. Backlash, ^ , 



21. Print consecutive number of this parameter run and all the 



parameters used: 



A 


J yi. 


V* 


e 


A 


/> 


J m 


Jl 


Jm + /* £ Jl 


fm 


fu 


r* 


K 


V Jn ?c 






£>» 


f 



22. Set program so always acts as combined system. 

25. Print consecutive number of this parameter run and selected 
parameters for linear syston.: 



2j e, 6) w =A Jl K 






2k. Solve the second order equation by Runge-Kutta-Gill numerical 

integration every 0.01 seconds problem time when zeta is greater 

than 0.1, otherwise every 0.004 seconds. Initial conditions 

• • • 

are 1*0, ©c= ©c = 0 and ©q - 1*0 or the pre- 

viously computed point. Next point is computed by four 
iterations of: 

~ ^ h (®R ~ ©c) ” & j ^ CO h ©c 

25. ’ Set EXIT from print routine. 

* 

26. Load time, an ^- P r i n ^ "buffer. 

27. Go to LIMITCYC to print only the overshoot. Otherwise print 
the phase trajectory points for each 0.1 seconds up to a maxi- 
mum of 12 cycles. 

m 

28. Store 0 for and © m 
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29. 

30. 




32. 



33. 



3^. 



35. 

36. 



37. 



38. 



39. 

40. 

41. 

42 . 

43. 

44. 
45- 
4 6 . 



Is problem time 0.1 sec.? 

Is problem time greater than 0.1 sec.? 

Print the number of point computed, time, ©^ and €> c . 

% 

EXIT (automatically set for desired Jump-out.) 

Store time, ©c , and for motor and load initial 

conditions. 

Is slope =? Ng? 



*>„+ ik. . 
, @C J L 

Is (5> negative? 

£ 

Store © c - A in 0 m 



- positive value? 



Solve 



©c= 



for & and ©£ 



by Runge-Kutta-Gill. 




for and ©*»by Runge-Kutta-Gill 

0 L ©m < o ? 

© L © ** ~“ A = negative? 

Set EXIT from print routine. 

Set EXIT from print routine. 

* 

Load time, ©e and 0c • 

Load 0^ and (0m • 

Print number of point computed and stored values. 
EXIT (automatically set for desired Jump-out). 
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k'J. Solve 

6 / 

©e ' 



*/ 

© 



J w +/>*Jl Y. Jn> 
Q - e(/g<3r* + ®c) 



“ e ) 



rv» 



OO / 



"Pm 

J *■« 



O/ 

©m + 



-y— (©«. " ©<) 

JlVI 



48. 

49. 

50. 

51. 

52. 
55. 

54. 

55. 

56. 

57. 

58. 

59. 

60. 

61. 

62. 
65. 

64. 

65. 

66 . 



Set EXIT from print routine. 

0 

Is ©c positive? ' 

Set up to count one phase trajectory when ©« is negative. 

16 9c > 1 - 0? 

Save 12 cycles of phase trajectory "been ccmplet ed? 

Arrange for no counting until 0 C is positive again. 

• • 

Is ^ — ©c. negative? 

Set up — ~ S>c) . 

% 0 

Is yOS^\ m greater than preselected £ ? 

Set up combined system. 

Xs ©c positive? 

Is negative? 

Is ) positive? 

Is 0^ — | negative? 

Is Oc-| negative? 

®c. y»©w= 01 
©e*/°©m=‘ 01 

Xs positive? 

Is 9 C ^^positive? 
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67. 



68 . 



6 9 . 



70. 



71. 



*72. 



75. 

74 . 



75. 

76. 



77* 



78. 



79. 



80. 



81 . 



82 . 



85- 



84. 



85- 



86 . 



87- 



Is Slope N l = Ng ? 

(l - Sc) C On 



O 

9c 



- Sf + 

* 

Modify FL0ATE5T to keep system 
vhen e c is negative. 

Is 0£ — 0^ positive? 

Store 0 L. “ A in @ ^ 




combined. 



=? positive value? 

Takes off modification 



Store ©C in e W\ . 

/* 

Is 0c - I negative? 

©C 

0c “ ©w = o? 

IS 0£ ~ ^»^P OSl - t:I - Ve ^ 

Is 0c -y/O 0** positive? 

Modify FLQATEST to keep system combined. Takes off modification 

c 

when Oc is positive. 

Have ll/2 minutes of real time elapsed since the last print-out? 
Is 0c < 1.0? 

Set up to print when in the fourth quadrant. 

Check if in the fourth quadrant. 

Is the new ©e greater than the previous 

G 

Store the new time, ©c and @ e . 

EXIT (automatically set for desired jump-out),. 

Print the last overshoot stored in UXETBUF (symbol 888 ) 

Stop real time clock. 

Print position in phase trajectory where it was stopped by the 
real time clock (symbol 999 ). 
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88 



. Return to beginning to start new problem if new parameters 
remain to be evaluated. Unconditional stop if all parameters 
have been evaluated. 

% 

89 . Has the required number of overshoots been evaluated prior to 
starting the averaging routine? 

90 . Add the <$£ just evaluated to the previous sum for averaging. 

91. Add one to the count prior to averaging. 

92. Clear time clock. 

o 

93. Print time, and ®c. • 

9^. Required number of total print-outs? 

95* Clear buffers. 

98 . Obtain the average of . 

97. Subtract 1.0. 

98 . Print the average of — 1.0 and how many overshoots were 
• utilized to obtain the average. 

99* Is, less than 1.005 radians? 

100. Clear buffers. 

101. Start real time clock. 

102. KX1T (automatically set to desired jump-out). 

103* Has 30 seconds of real time elapsed since the last print-out? 
10lf. Print the position in the phase trajectory where it was stopped 
by the real time clock and < 3 : being less than 1.005 radians 
(symbol 1005 ) • 
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Jump Switches: 



- For LBUTCYC , allows the maximum overshoot print-out only. 

With J down, a phase trajectory is printed for a maximum 

« 

of 12 cycles. 

J^ - Set the number of print-outs depending on zeta. Change 
the number of overshoot values used to obtain the average 
value. 

J - Use the real time clock for an automatic recycle to a 

3 

new parameter when have: 

a. ll/2 minutes maxi mm after the last print-out. 

b. 30 seconds maximum and is 1.005 radians or less. 



Stop Switches: . 

-S^ - Stops at the end of this parameter run or stops at the end 
of the manual print-out. 
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APPENDIX B COMPUTER PROGRAM 











ORG 




00007 


00007 


74 


0 


00070 

40701 


fault exf 


0 


00070 




75 


0 


SLJ 


0 


FAULPRIN 










ORG 




40000 


40000 


74 


0 


00100 


start exf 


0 


00100 


40001 


11 


0 

1 


40004 

40623 


LAZYZETA L^ 


0 

1 






20 


0 


40425 


STA 


0 


ZETA 


40002 


50 


1 


00000 


ENI 


1 


0 


75 


0 


40000 


SLJ 


0 


START 


40003 


75 

50 


0 

0 


SoSoo 


MANSTART SLJ 
ENI 


0 

0 


MANUAL 

0 


40004 


12 


0 


40565 


LDA 


0 


newstop 




20 


0 


40563 


STA 


0 


STOPNOW 


40005 


12 


0 


40425 


NEWZETA LDA 


0 


ZETA 


22 


0 


40564 


A JP 


0 


STOP I NDX 


40006 


32 


0 


40665 


FMU 


0 


TWO 


40007 


32 


0 


40426 


FMU 


0 


OMEGAN 


20 


0 


40404 


STA 


0 


A 


12 


0 


40624 


LDA 


0 


ZETAINDX+1 


40010 


31 


0 


40425 


FSB 


0 


ZETA 


22 


0 


40433 


A JP 


0 


QUIKCOMP 


40011 


22 


2 


40433 


AJP 


2 


QUIKCOMP 




12 


0 


40402 


LDA 


0 


P0INT01 


40012 


20 


0 


40154 


STA 


0 


TABLES+1 


40013 


20 

20 


0 

0 




STA 

STA 


0 

0 


TABLEL+1 

TABLEM+1 


12 


0 


40661 


LDA 


0 


COUNTlO 


40014 


20 


0 


40400 


STA 


0 


tenthsec 


50 


0 


00000 


ENI 


0 


0 


40015 


75 


2 


40437 


ZZ SLJ 


2 


BIGZETA 


50 


0 


00000 


ENI 


0 


0 


40016 


12 


0 


40652 


LDA 


0 


STOP20 


40017 


20 

12 


8 


48 U 2 


STA 

LDA 


0 

0 


STOPZETA 
COUNT 12 




2° 


0 


40617 

40672 


STA 


0 


COUNTZET 


40020 


12 


0 


LDA 


0 


EIGHT 


20 


0 


40620 


STA 


0 


DIVIDE 


40021 


10 


0 


00010 


ENA 


0 


10 


61 


0 


40562 


SAL 


0 


STOPRINT 


40022 


75 


4 


71000 

00000 


GOAROUND SLJ 


4 


DECOF 


50 


0 


ENI 


0 


0 


40023 


02 


0 


40427 


02 


0 


FORZEROS 


04 


0 


00000 


04 


0 


0 


40024 


12 


0 


40664 


LDA 


0 


ONE 


20 


0 


40156 


STA 


0 


UDOT 


40025 


10 


0 


00000 


ENA 


0 


0 


20 


0 


00000 


STA 


0 


0 


40026 


20 


0 


40162 


STA 


0 


THETA 


20 


0 


40161 


STA 


0 


THETADOT 


40027 


20 


0 


40157 


STA 


0 


U 


20 


0 


40155 


STA 


0 


T 


40030 


20 


0 


40377 


STA 


0 


INDEX 


20 


0 


40615 


STA 


0 


INDEXAVE 


40031 


20 


0 


40351- 


STA 


0 


PRI NTBUF+3 




20. 


0 


40352 


STA 


0 


PRINTBUF+4 


40032 


20 


0 


40353 


STA 


0 


PRI NT8UF+5 


20 


0 


40614 


STA 


0 


LIMINDEX 


40033 


20 


0 


-40607 


STA 


0 


LIM I TBUF+2 




20 


0 


40610 


STA 


0 


LIMITBUF+3 


40034 


20 


0 


40611 


STA 


0 


LIMITBUF+4 

LIMITBUF+5 


20 


0 


40612 


STA 


0 


40035 


20 


0 


40613 


STA 


0 


LIMITBUF+6 


" 


12 


0 


40621 


LDA 


0 


PRINTSET 


40036 


20 


0 


40326 


STA 


0 


OK2PRINT+1 




20 


0 


40335 


STA 


0 


BOUNPRIN+4 


40037 


12 


0 


40622 


LDA 


0 


NEWLIMPT 




20 


0 


40544 


STA 


0 


LIMPRINT+1 


40040 


12 


0 


40253 


LDA 


0 


FLOTDATA 




20 


0 


40071 


STA 


0 


FLOATEST 


40041 


12 


0 


40254 


■ LDA 


0 


FLOTDATA+1 


20 


0 


40072 


STA 


0 


FLOATEST+1 


40042 


75 


0 


40452 


NEWINPUT SLJ 


0 


I NDXVALU 


50 


0 


00000 


ENI 


0 


00000 
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40043 

40044 

40045 

40046 

40047 

40050 

40051 

40052 

40053 

40054 

40055 

40056 

40057 

40060 

40061 

40062 

40063 

40064 

40065 

40066 

40067 

40070 

40071 

40072 

40073 

40074 

40075 

40076 

40077 

40100 

40101 

40102 

40103 

*0104 

40105 

40106 

40107 



ll 


8 


*8148 






8 


Efhl&R 


12 


0 


40513 




LDA 


0 


NOS I N< 


20 

il 


0 

8 


8?8Xl 

00000 


. 


STA 

fM 


0 

4 

0 


floatest 

DECOF 

0 




0 


40404 


• 


01 


0 


A 


06 


0 


00001 




06 


0 




72 

75 


0 

4 


40046 

71000 


• 


RAO 

SLJ 


0 

4 


deJof 


01 


0 


40412 




01 


0 


MOINERT 


06 


0 


00000 




06 


0 


0 


72 


0 


40510 




RAO 


0 


LINEAR+10 


75 


4 


71000 




SLJ 


4 


DECOF 


01 


0 


40420 




01 


0 


KMOTCONS 


06 


0 


00000 




06 


0 


0 


10 


0 


00000 


combined 


ENA 


0 


0 


20 


0 


40160 




STA 


0 


OX 


20 


0 


40163 




STA- 


0 


QY 


75 


•4 


60200 




SLJ 


4 


RUNGE 


00 


0 


40153 




0 


0 


TABLES 


00 


0 


40164 




0 


0 


DERIVES 


75 


4 


60201 




SLJ 


4 


RUNGE+1 


50 


0 


00000 




ENI 


0 


0 


75 


0 


40310 




SLJ 


0 


COMBPRIN 


50 


0 


00000 




ENI 


0 


0 


12 


0 


40155 




LDA 


0 


T 


20 


0 


40174 




STA 


0 


TL 


20 


0 


40210 




STA 


0 


TM 


1 2 


0 


40156 




LDA 


0 


UDOT 


20 


0 


40175 




STA 


0 


VDOT 


33 


0 


40411 ' 




FDV 


0 


RHO 


20 


0 


40211 




STA 


0 


WDOT 


12 


0 


40157 




LDA 


0 


U 


20 


0 


40176 




STA 


0 


V 


20 


0 


40200 




STA 


0 


THETADL 


33 


0 


40411 




FDV 


0 


RHO 
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